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We study numerically the behavior of the Biham-Middleton-Levine traffic model in three dimen- 
sions. Our extensive numerical simulations show that the phase diagram for this model in three 
dimensions is markedly different from that in one and two dimensions. In addition to the full speed 
moving as well as the completely jamming phases, whose respective average asymptotic car speeds 
{v) equal one and zero, we observe an extensive region of car densities p with a low but non-zero 
average asymptotic car speed. The transition from this extensive low average asymptotic car speed 
region to the completely jamming region is at least second order. We argue that this low speed 
region is a result of the formation of a spatially-limited-extended percolating cluster. Thus, this low 
speed phase is present in n > 3 dimensional Biham-Middleton-Levine model as well. 
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I. INTRODUCTION 

With the ever increasing computational power, simu- 
lating trafhc in the microscopic level by means of cellular 
automaton becomes a real possibility. One of the sim- 
plest model for city traffic of this kind is the so-called 
Biham-Middleton-Levine (BML) traffic model [Q . 

The one-dimensional BML model is simply the elemen- 
tary binary CA rule 184 operating on a one-dimensional 
lattice with periodic boundary condition. The asymp- 
totic car speed (v) in this one-dimensional model is ex- 
actly known, and is given by 



{v) = 



1 



if p< 1/2, 
if 1/2 < p < 1, 



(1) 



where p is the car density in the system l^]. In other 
words, in the one dimensional BML model, traffic jam 
occurs only when the car density p is equal to 1 = pi^' 
and all cars move in full speed whenever p < 1/2. 

The two-dimensional BML model considers the mo- 
tions of north- and east-bound cars in a two-dimensional 
square lattice with periodic boundary conditions in both 
the north-south and east-west directions. (We shall give 
the exact rules for the two-dimensional BML model in 
Section ||.) Although we lack an exact analytical expres- 
sion for the average asymptotic car speed (v) as a func- 
tion of car density p in the two-dimensional BML model, 
extensive numerical simulations ||l| as well as mean field 
theory studies j^] have been carried out. Their results 
strongly suggest a fluctuation-induced first order phase 
transition in (v). Moreover, the average asymptotic car 
speed is likely to follow 



{v) = 



for < p < pi^ 

f ^ (2) 

for p > p^ % 



(2) 

where the critical density pc is numerically found to be 
about 0.31 [Q and is analytically proven to be less than 
1/2 In addition, Tadaki and Kikuchi found a more 
subtle phase transition related to the final jamming pat- 
tern. Their numerical study showed that jamming pat- 
terns for car density p less than about 0.52 are very well 
self-organized. On the other hand, when p is greater than 
0.52, the jamming patterns are random 

Extension of the BML model to higher dimensions can 
be regarded as a highly simplified model for computer 
network communication in a hypercube. And from the 
physics point of view, it is natural to investigate the phase 
diagram as well as the upper critical dimension of the 
BML model in higher dimensions. As a pioneer study, we 
report the result of an extensive numerical study of the 
BML model in three dimensions in this paper. We find 
that the three dimensional model has a richer phase dia- 
gram than that in one and two dimensions. In addition 
to the fiuctuation-induced first order phase transition in 
(v), we also observe a low (but non-zero) speed phase. 

To begin, we first introduce the higher dimensional 
generalization of the BML model in section ||. Then, we 
report our simulation results in section |l| and present 
our analysis of results in section IV. Finally, we draw 
our conclusions in section 



(2) 



II. THE BML MODEL 



Let us introduce the modified BML model in n- 
dimensions below. Consider an n-dimensional NiX N2X 
• • • X Nn square lattice with periodic boundary condi- 
tions. Each lattice site will either contain no car (that is. 
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an empty site) or contain exactly one car moving in the 
direction. We denote pi the density of cars moving along 
e-i. (That is, pi equals the number of cars moving along 
the Ci direction divided by the total number of cars in 
the system.) We denote the total car density of the sys- 
tem hy p = "Y^ Pi and we define the car density vector by 
p = p2, ■ • ■ , Pn)- Initially, cars are placed randomly 
and independently onto the n-dimensional square lattice 
according to a pre-determined car density vector p. 

The dynamics of the cars are governed by the following 
rules. Each ei-moving car advances one site along the ei 
direction provided that no car blocks its way. Otherwise, 
that ei-moving car stays in its present location. Parallel 
update is taken for all ei-moving cars. After this, each 
e2-moving car advances one site along the 62 direction 
if no car blocks its way. Otherwise, that e2-moving car 
stays in its present location. Again, parallel updating 
is used. This process goes on until each e„-moving car 
is given a chance to move. This marks the end of one 
timestep and the above car moving process is repeated 
over and over again. 

At each timestep, the car speed is defined as the ratio 
of number of cars moved to the total number of cars in the 
lattice system. And the average asymptotic car speed {v) 
is defined as the car speed averaged over both the cycle 
time and initial car configurations. Since we are inter- 
ested in the behavior of the system in thermodynamic 
limit, so we only consider the limit when Ni, N2, ■ ■ ■ , Nn 
all tend to infinity while keeping the aspect ratio between 
each side fixed. 

We define the n-dimensional BML traffic model to be 
the one with aspect ratio between each of the n sides be- 
ing fixed to one. That is to say, A'^i = ^"2 = • • • = Nn- 
Also, an n-dimensional BML traffic model is called ho- 
mogeneous if and only if p,; = pj for all i,j In 
this paper, we concentrate on the homogeneous three- 
dimensional BML traffic model. So for simplicity, we 
shall simply call it the three-dimensional BML model 
whenever confusion is not possible. From this definition, 
it is clear that the average asymptotic car speed {v) is 
a function of p only and its value lies between zero and 
one. 



III. SIMULATION RESULTS OF THE 
THREE-DIMENSIONAL BML MODEL 

Our simulation is performed on a variety of machines 
including clusters of Sun Sparc and Dec Alpha worksta- 
tions, various Pentium-based PCs and Power PCs as well 
as the SP2 supercomputer. The estimated total CPU 
time is about 300 mips years. Even so, owing to our CPU 
time limitations, we can only systematically simulate up 
to a lattice size of 100 x 100 x 100. Nonetheless, we have 
also simulated for the cases of very small and very large 
car densities up to a lattice size of 1000 x 1000 x 1000 
before finally drawing our conclusions. 
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FIG. 1. The car density p vs. average asymptotic car speed 
{«) for the three-dimensional BML model on a 100 x 100 x 100 
lattice. In the two inserted blowups, one clearly observes 
the drop of (,v) from 1 to NI{N -|- 1) at p = 0.01 as well 
as a very small bump just below the first critical density 
Pci' = 0.18 ±0.01. Besides, one also observes that (u) reaches 

(3) 

zero at the critical car density p^^ = 0.32 ± 0.02. 

Fig. |l| shows the p vs. (w) curve for the BML model 
in a 100 x 100 x 100 lattice. Each data point in the fig- 
ure represents the average asymptotic car speed over an 
ensemble of random initial configurations. For p < 0.1 
as well as p > 0.22, the value of (v) is obtained by 
averaging over 1000 initial configurations. In contrast, 
(v) for 0.1 < p < 0.22 is obtained by averaging over 
only 50 random initial configuration because the long 
relaxation time prevents us from obtaining more sam- 
ples. Fig. tells us that (v) — 1 when the car density 
p < 0.005 « l/2iV. We call this region the "fuU speed 
phase" . Moreover, recurrent states in the car density re- 
gion are cycles of period 100 = N. (The dependence of 
various parameters on N here and hereafter are based 
on our simulation results in various lattice sizes up to 
1000 x 1000 X 1000, including various odd, even and prime 
values of N.) As p increases to about 0.01 ~ 1/iV, {v) 
begins to drop. The recurrent states form cycles with 
periods several ten times the linear lattice size N. As 
the car density p reaches about 0.02 « 2/A^, {v) drops 
to the value N/{N + 1) ^ 100/101 « 0.99099 and stays 
constant until p reaches about 0.10. In this car density re- 
gion, recurrent states form cycles of period 101 = + 1. 
In other words, in the recurrent state, each car in the 
system will be blocked exactly once in each cycle. For 
the 100 X 100 X 100 lattice, (v) is slightly greater than 
N/{N + 1) as 0.10 < p < 0.17. A similar but much 
smaller bump in {v) is also observed in the 200 x 200 x 200 
lattice. Hence, we conclude that the bump is due to fi- 
nite size effect of the lattice. Typical recurrent states 
in this car density region are cycles of period « lOOA^. 
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Moreover, the typical relaxation time for a random con- 
figuration in this range of car density appears to scale 
exponentially with N. In fact, the long relaxation time 
forbids us from performing systematic simulation with 
lattice size greater than 100 x 100 x 100. Clearly, the 
exponentially long relaxation time signals a critical slow- 
down. 

As the car density reaches about — 0.18 ±0.01, (v) 
drops abruptly to about 0.015. Since the asymptotic car 
speeds in all our simulation data are either greater than 
or equal to N/{N + 1) or less than 0.03, we strongly be- 
lieve that the observed sudden drop in {v) is a result of a 
first order phase transition. Interestingly, the periods of 
the recurrent configurations of all these low but non-zero 
speed states are equal to N. When we further increase 
the car density p, (v) gradually decreases until it finally 

reaches zero at p^^' = 0.32 ± 0.02. 

In summary, our simulation tells us that for a finite 
N X N X N lattice, the system exhibits a non-trivial 
"high speed region" with (v) = N/{N + 1) as well as 
a non-trivial "low speed region" with (v) < 0.03. Thus, 
in the thermodynamic limit, the three-dimensional BML 
model has a full speed phase, a low speed phase and a 
completely jammed phase. (Moreover, just like the two- 
dimensional case, the completely jammed phase may fur- 
ther be divided into the self-organized jamming and the 
random jamming regions.) The transition from the full 
speed to the low speed phase is first order in nature and 
the transition from the low speed phase to the completely 
jammed phase is smooth. That is to say, this transition 
is at least second order. 




FIG. 2. A typical high speed phase recurrent configuration 
of period six in a 5 x 5 x 5 lattice. Readers may verify that each 
car in this configuration is blocked exactly once per period. 



IV. ANALYSIS OF OUR SIMULATION RESULTS 



A. The Full Speed Phase 



In addition to the systematic trend that the critical 
car density from the high to low speed phase strictly de- 
creases with spatial dimension, we observe an interesting 
feature in the high speed phase of the three-dimensional 
BML model. Unlike the one- and two-dimensional mod- 
els, the recurrent configurations for any finite N x N x N 
lattice in the high speed region with p^l/N form cycles 
of period + 1. A typical high speed recurrent config- 
uration in a 5 X 5 X 5 lattice is shown in Fig. |^ as an 
illustration. Readers may verify that cars in these high 
speed configurations will be blocked once per cycle period 
0. Unfortunately, we do not have a good explanation 
why this {v) — N/{N + 1) recurrent state is preferred 
over the {v) ~ 1 recurrent state in three dimensions. 



B. The Transition To Low Speed Phase 

Since we cannot find any intermediate speed asymp- 
totic configurations in our simulation, the transition from 
the high to low speed phase is likely to be first order. To 
further investigate to the nature of this transition, we 
drive the system by slowly adding cars to or removing 
cars from the system. That is to say, starting from p — 0, 
we increase the car density by a fixed small amount Ap 
by randomly introducing cars to the empty sites in the 
system. And then, we evolve the system until it relaxes to 
a recurrent state. We repeat the process until p reaches 
one. After this, we decrease the car density of the system 
by Ap by randomly removing cars from the system. And 
then, we evolve the system until it reaches a recurrent 
state. We repeat this process until p becomes zero. The 
p vs. (v) graph obtained in this way on a 100 x 100 x 100 
lattice with Ap = 0.001 is shown in Fig. ^ Clearly, as 
we slowly increase the car density p, transition to the low 
speed phase occurs at car density around 0.22, which is 

(3) 

slightly higher than the critical car density pc/ « 0.18. 
More dramatically, as we slowly decrease the car density, 
transition to the high speed phase occurs at car density 
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around 0.07, which is much smaller than the critical car 

(3) 

density /Oci • The observed hysteresis loop confirms the 
hypothesis that this is a first order phase transition. And 
since the only nonlinearity in the model comes from the 
exclusion volume effect, we conclude that the phase tran- 
sition is fluctuation induced. 
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FIG. 3. A typical evolution of car density p vs. asymptotic 
car speed (v) in a 100 x 100 x 100 lattice when we drive the 
system slowly as discussed in the text. The solid and dash 
curves represent the evolution of the system when cars are 
slowly added to and removed from the system, respectively. 

C. The Low Speed Phase And The Completely 
Jamming Phase 

Unlike its one- and two-dimensional counterparts, the 
three-dimensional BML model has a low speed phase 
with < (v) ^ 0.03. Similar to the completely jam- 
ming configurations, we find that the recurrent config- 
urations in the low speed phase contain (directed) per- 
colating clusters of cars. But unlike the completely jam- 
ming configurations, we found a small number of residual 
freely moving cars in the low speed phase. Hence, the pe- 
riod of these recurrent states equals the linear system size 
TV. And since most cars are already jammed by colliding 
into the percolating cluster, the average asymptotic car 
speed is low. A typical low speed recurrent configuration 
in a 5 X 5 X 5 lattice is shown in Fig. ^ as an illustration 

Recall that a percolating backbone is essentially an 
one-dimensional object. Therefore, if the background lat- 
tice is one- or two-dimensional, all other moving cars will 
eventually merge into the percolating backbone leading 
to a completely jamming configuration. The situation 
is completely different when the background lattice is at 
least three-dimensional. In this case, since both the tra- 
jectories of moving cars and the percolating backbone are 
essentially one-dimensional objects, trajectories of some 



moving cars may not intersect with the percolating clus- 
ter. Thus, if the car density is small enough, some resid- 
ual freely moving cars may present in a recurrent config- 
uration giving rise to the observed low speed phase. 

As the car density gradually increases in the low speed 
phase, the size of the percolating cluster in the recur- 
rent configuration increases. It becomes more and more 
difficult for the system to accommodate residual freely 
moving cars. Hence, (v) gradually decreases until it even- 
tually reaches zero. The transition from the low speed 
phase to the completely jamming phase is, therefore, 
smooth. 




FIG. 4. A typical low speed recurrent configuration in 
a 5 X 5 X 5 lattice making up of five full speed cars with 
{v) — 1/14. We denote the five full speed cars by bold ar- 
rows. 



V. CONCLUSIONS AND OUTLOOK 

In summary, we study the phase diagram of the three- 
dimensional BML model. Similar to the two-dimensional 
model, a fiuctuation-induced first order phase transition 
in asymptotic average car density (v) is observed at car 
density pc/ = 0.18 ± 0.01. We also discover a new low 
speed phase which is absent in one- and two-dimensional 
models. We argue that the existence of this low speed 
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phase is geometrical in nature, and hence this phase 
will exist in higher dimensional BML models as well. 
It is instructive to numerically verify our claims in the 
four-dimensional model. Unfortunately, the amount of 
computation involved will probably be too high for us 
at this moment. Finally, our simulation suggests that 
the transition from the low speed phase to the com- 
pletely jamming phase is smooth and occurs at car den- 
sity pc^ = 0.32 ± 0.02. 

A number of open questions remain. For instance, we 
do not understand why the {v) = N/{N + 1) high speed 
states are preferred over the (v) = 1 full spc;ed states 
in the three-dimensional model. And it is meaningful to 
investigate if this behavior persists in higher dimensions. 
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